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ABSTRACT 

Here we present a self-consistent stationary solution for spherically symmetric winds driven 
by massive star clusters under the impact of radiative cooling. We demonstrate that cooling 
may modify drastically the distribution of temperature if the rate of injected energy approaches 
a critical value. We also prove that the stationary wind solution does not exist whenever the 
energy radiated away at the star cluster center exceeds ~ 30% of the energy deposition rate. 
Finally we thoroughly discuss the expected appearance of super-star cluster winds in the X- 
ray and visible line regimes. The three solutions here found: the quasi-adiabatic, the strongly 
radiative wind and the inhibited stationary solution, are then compared to the winds from Arches 
cluster, NGC 4303 central cluster and to the supernebula in NGC 5253. 
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In the stationary solution for spherically sym- 
metric winds (Chevalier and Clegg 1985; hereafter 
referred to as CC85) as well as in the former ap- 
proach of Holzer and Axford (1970) and in the 
more recent numerical calculations of Canto et al. 
(2000) and Raga et al. (2001) the flow has been 
assumed to be adiabatic and thus predicts a very 
extended X-ray envelope around the sources. The 
impact of cooling on the stationary wind solution, 
was discussed by Silich et al. (2003, hereafter re- 
ferred to as Paper I) for winds driven by powerful 
and compact stellar clusters, and by Wang (1995) 
for gas outflows from galaxies. Winds driven by 
compact star clusters establish a temperature dis- 
tribution radically different from that predicted by 



the adiabatic solution, bringing the X-ray emitting 
boundary much closer to the star cluster surface. 
However, in none of the above studies, the effects 
of radiative cooling within the star forming volume 
itself were taken into consideration. 

Here we present a self-consistent semi-analytical 
model of stationary winds driven by massive stel- 
lar clusters taking full account of radiative cooling 
(see sections 2 and 3). We first discuss how to 
find proper wind central values and then use them 
to integrate numerically the basic equations. We 
also indicate the threshold value of the energy de- 
position rate above which a stationary solution is 
inhibited. In sections 4 and 5, the three regimes 
found when radiative cooling is considered: the 
quasi-adiabatic, the strongly radiative wind and 
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the inhibited stationary wind, are then compared 
to well observed examples. Our conclusions are 
given in section 6. 

2. The adiabatic solution 

Following CC85, Canto et al. (2000) and Raga 
et al. (2001), we assume that within a star clus- 
ter, within the volume of radius R SCl the matter 
ejected by stellar winds and supernova explosions 
is fully thermalized via random interactions. This 
generates the large central overpressure that con- 
tinuously accelerates the ejected gas and eventu- 
ally blows it out of the star cluster volume. There 
arc three star cluster parameters which together 
define the hydrodynamical properties of the resul- 
tant wind outflow (or the run of density, tempera- 
ture and expansion velocity, which asymptotically 
approach p w ~ r~ 2 , T w ~ r~ 4 / 3 u w w V^a)- 
The three parameters are: the total energy (E sc ) 
and mass (M sc ) deposition rates and the actual 
size of the volume that encloses the star clus- 
ter (i? sc ). The total mass and energy deposi- 
tion rates also define the wind terminal velocity 

Vooa = \J2E SC /M SC . 

In the adiabatic case there is an analytic so- 
lution and thus one can derive the wind central 
density, pressure and temperature (see Canto et 
al., 2000) if the above parameters are known: 
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and q m are the energy and mass deposition 
rates per unit volume (q e = 3E ' sc / 'Air R sc \ q m = 
3M sc /47ri?g C ), 7 is the ratio of specific heats, n 
is the mean mass per particle and k is the Boltz- 
mann constant. Using these initial values one can 
solve the stationary wind equations numerically 
and reproduce the analytic solution throughout 
the space volume. The relevant equations are: 
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In equations (4-6) r is the spherical radius and Q 
is the cooling rate, assumed equal to zero in CC85, 
Canto et al. (2000) and Raga et al. (2001). 

Within the central volume, temperature and 
density present almost homogeneous values, where- 
as the expansion velocity grows almost linearly 
from km s _1 at the center, to the sound speed 
at the cluster radius, r = R sc . There is then a 
rapid evolution as matter streams away from the 
star cluster. The flow accelerates rapidly when 
approaching the sonic point and the wind tem- 
perature and density begin to deviate from their 
central quasi-homogeneous distributions. At large 
radius the resultant wind parameters rapidly ap- 
proach their asymptotic values (see CC85 and 
paper I). 

3. The radiative solution 

Due to the highly nonlinear character of the 
cooling function, the analytic approach is not valid 
in the general case that includes radiative cooling, 
and thus one needs to perform a numerical inte- 
gration. However in such a case, the star clus- 
ter parameters (E sc , M sc and R sc ) do not define 
the wind central temperature and density and the 
problem arises: how to solve equations (4-6) if nei- 
ther the initial nor the boundary conditions are 
known? 

To solve the problem we re-write equations (4- 
6) and obtain within the star cluster radius (r < 

Rsc) 

du w = 1 (l-l)(qe-Q) + q m (^l-lcl) 
dr p w c 2 - ul, 
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and for r > R sc 
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In equations (7-12) P w is the wind thermal pres- 
sure, c s — (7P/ p) 1 ! 2 is the sound speed, Q — n^A, 
n w is the wind atomic number density and A(Z, T) 
is the cooling function (a function of metallicity 
and temperature, see Raymond et al. 1976). 

It is easy to prove that the derivative of the ex- 
pansion velocity is positive throughout the space 
volume, only if the sonic point is at the star cluster 
surface. Indeed, if the expansion velocity becomes 
supersonic at r < R sc , the right-hand side of equa- 
tion (7) and du w /dt become negative inside the 
star cluster volume. On the other hand, if the ex- 
pansion velocity remains subsonic at r > R sc , the 
right-hand side of equation (10) and the derivative 
of the expansion velocity, become negative outside 
the star forming region. 

The above implies that a stationary wind so- 
lution, which assumes a continuous gas acceler- 
ation, exists only if the outflow crosses the star 
cluster surface at the local sound speed (u = c s at 
r = R ac ). This conclusion, illustrated in Figure 1 
doesn't depend on the wind thermodynamic prop- 
erties. It is valid both for the adiabatic (CC85 and 
Canto et al. 2000) and for the radiative solution. 

There are three possible types of integral curve 
(Figure 1) corresponding to different possible po- 
sitions of the sonic point with respect to the star 
cluster surface. 1) The stationary wind solution 
(solid line, R son ic = Rsc)- In this case the ther- 
mal pressure decreases continually outside the star 
cluster surface, approaching a negligible value at 
large radii. 2) The breeze solution (dashed line, 
Rsonic > Rsc)- In this case the central temper- 
ature T c is smaller than in the first case, which 
shifts the sonic point outside the cluster. This 
branch of solutions requires of a finite confining 
pressure and leads to zero velocity at infinity, in 
agreement with Parker's (1958) conclusion. 3) The 
unphysical double valued solution (dotted line, 
Rsonic < Rsc)- In this case T c is larger than in the 
stationary wind case. The thermal pressure goes 
to zero when the expansion velocity approaches 
the wind terminal speed. 

The appropriate solution is selected by the cen- 
tral conditions. In order to obtain a stationary 




Fig. 1. — Three possible types of integral curves. 
1) The stationary wind solution (solid line), 
Rsonic = Rsc- 2) The breeze solution (dashed 
line), Rsonic > Rsc- 3) Unphysical double val- 
ued solution (dotted line), R SO nic < Rsc- For the 
three examples we adopted E sc = 10 41 erg s _1 , 
R sc = 15 pc, Vw = {2q e /q m ) 1 ' 2 = 1000 km s" 1 . 



free wind solution, one has to find the wind cen- 
tral density and central temperature which accom- 
modate the sonic point at the star cluster surface 
(u = c s at r — R sc )- This is the key point that 
allows for the definition of the central (r = 0) wind 
parameters and for the numerical solution of equa- 
tions (7-12). 

The wind central temperature T c and central 
atomic number density n c are not independent in 
the radiative case. They are related by the equa- 
tion 
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This results from comparing the derivative of the 
expansion velocity at the star cluster center, using 
equation (9), with equation (7). Note that in the 
absence of radiative cooling, when 



Q = n 2 c A{T c ) = q e - 



7- 



< = 0, 



(14) 



equation (13) is transformed into the adiabatic re- 
lation (3). Thus, the wind parameters at the star 
cluster center can be found by iteration of the cen- 
tral temperature until the sonic point takes its 
proper position at the selected star cluster surface. 
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As in the adiabatic case p c and T c are indepen- 
dent, one can always find the central pressure that 
accommodates the sonic point at the star cluster 
surface (see equations 1 and 3): 



67 



BP r . 



7 - 1 V 2 <7e<7r, 



(15) 



However this is not the case if radiative cooling 
is taken into account. In this case the central 
density and the central pressure go to zero when 
the central temperature approaches its maximum 
value T max = ^ (see equation 13). P c in- 

creases for smaller values of the central tempera- 
ture. However, it cannot exceed a maximum value 
(see equation 13 and Figure 2a), bound by the gas 
radiative cooling. At this critical stage the fraction 
of energy radiated away at the cluster center per 
unit time, 5 = (q e — n 2 A(T c ))/q e , reaches 30% 
of the injected energy (see Figure 2c). If the cen- 
tral temperature becomes even smaller, there is no 
density enhancement able to compensate the fall 
in pressure promoted by radiative cooling. Con- 
sequently, the central pressure cannot promote an 
effective outward acceleration. Therefore in the 
radiative case, the sonic radius (R SO nic) cannot be 
arbitrarily large and has a maximum value for any 
given set of star cluster parameters. 

Figure 2b shows how R SO nic depends on the cen- 
tral temperature, for particular values of q e and 
q m . For the largest central temperature, the resul- 
tant central pressure acquires its lowest value (see 
Figure 2a) and the sonic point lies very close to the 
star cluster center. Smaller values of T c lead to a 
larger P c (Figure 2a) and consequently to larger 
Rsonic values. Rsonic reaches its maximum pos- 
sible value when the central pressure approaches 
also its maximum (compare figures 2a and 2b) as 
expected in the quasi-adiabatic regime (see equa- 
tion 15). 

In the radiative stationary solution, although 
the sonic point may approach its maximum pos- 
sible value, radiative losses of energy would rep- 
resent only a moderate fraction of the energy in- 
put rate. In these cases, cooling would drastically 
modify the wind temperature distribution outside 
the star forming region (see Paper I and discussion 
below). When the rate of energy radiated away 
at the star cluster center exceeds <~ 30% of the 
energy deposition rate, the stationary solution is 
inhibited. This conclusion is stressed in Figure 3, 
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Fig. 2. — The impact of the central temperature 
on the outflow, a) The outflow central pressure; b) 
The position of the sonic point; c) Comparison of 
the deposited and the radiated energies at the star 
cluster center. In all cases the energy and the mass 
deposition rates per unit volume are q e = 2.4 x 
10~ 18 erg cm~ 3 s _1 and q m = 4.8 x 10~ 34 g cm~ 3 
s~\ respectively. (2q e /q m ) 1 / 2 = V^.a = 1000 
km s _1 , and the assumed ejected gas metallicity 
is solar. 
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Fig. 3. — The impact of radiative cooling. The 
threshold energy input rate above which the sta- 
tionary wind solution is fully inhibited, as function 
of the star cluster radius. The solid line repre- 
sents the threshold energy for star clusters with 
(2q e /q rn ) 1 / 2 = V xA = 1000 km s" 1 . The dotted 
and the dashed lines mark the threshold energies 
for star clusters with (2q e /q m ) 1 ^ 2 = 1500 km s _1 
and 500 km s _1 , respectively. 

which displays the critical energy deposition rate 
for different values of q e /q m = V^ oA /2. Moving 
from right to left along the horizontal line is equiv- 
alent to considering progressively more compact 
clusters, all with the same energy and mass depo- 
sition rates (E sc w 4.4 x 10 41 erg s _1 ;M sc w 1.4 
M Q yr _1 ). For large star clusters the maximum 
allowed sonic radius R. SO nic exceeds the star clus- 
ter radius R sc , however one can accommodate the 
sonic point at the star cluster surface once a proper 
central temperature is selected and obtain a sta- 
tionary wind solution. However, if the consid- 
ered star cluster is smaller than the critical value 
(~ 12 pc for the example shown in Figure 3), the 
maximum allowed sonic point radius moves inside 
the star cluster and the stationary wind solution 
vanishes. 

The same is true if one moves along the verti- 
cal line in Figure 3, from low to high energy input 
rates. In this case one is selecting progressively 
more energetic star clusters within the same vol- 
ume, until the sonic point ends up inside the star 
cluster (in our example at L cr u ~ 4.4 x 10 41 erg 
s _1 ) and the stationary wind solution vanishes. 

Once the proper initial conditions are selected, 



one can solve the main equations (7-12) numeri- 
cally and obtain the wind temperature and den- 
sity distributions. We have compared for example 
our results with Stevens & Hartwell (2003) stan- 
dard model (R sc = 1 pc, M sc = 1CT 4 Mgyr" 1 , 
VocA = 2000 km s _1 ). In this case the station- 
ary wind evolves in the quasi-adiabatic regime and 
we found an excellent agreement with Stevens & 
Hartwell central values and X-ray luminosity. Our 
model predicts T c = 5.9 x 10 7 K, n c = 0.65 cm~ 3 
and the X-ray flux between 0.3 and 8.0 keV from 
the central lpc volume L x = 5.2 x 10 32 erg s _1 . 

4. The expected appearance of stationary 
radiative winds 

Free winds present a four zone structure (Silich 
et al. 2003): a star cluster region filled with a hot 
X-ray plasma, an adjacent X-ray halo with a de- 
creasing temperature distribution, the line cooling 
zone and a region of rccombined gas, exposed to 
the UV and soft X-ray radiation from the inner 
zones and to the UV photons emitted by the star 
cluster itself. 
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Fig. 4. — The impact of cooling in the extended X- 
ray zone of stationary winds. Temperature profiles 
for progressively larger energy deposition rates. 
Solid, dotted and dashed lines represent winds 
with 0.5, 2 and 3 xlO 41 erg s _1 , respectively. 
R sc = 10 pc in all three cases here considered. 

Figure 4 presents the free wind temperature dis- 
tribution for three SSCs, all with the same R sc = 
10 pc radius and the same ratio (2g e /g m ) 1 / 2 = 
^00,4 = 1000 km s _1 , but different energy and 
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mass deposition rates. The lowest energy case 
(solid line) lies well into the adiabatic regime. In 
the other two cases however, the radiative cooling 
clearly modifies the internal wind structure bring- 
ing the boundary of the X-ray zone and the pho- 
toionized envelope closer to the star cluster sur- 
face. In the most energetic case shown in Fig- 
ure 4, that with 3 x 10 41 erg s _1 star cluster 
(dashed line), the outer boundary of the X-ray 
zone (Tx-ray ~ 5 x 10 5 K) is about a factor of 
1.5 smaller than in the adiabatic case. Further- 
more, the dimension (R4) at which the gas attains 
a temperature ~ 10 4 K lies about 10 times closer 
to the star cluster center than in the adiabatic 
case. Consequently, the maximum density of the 
emission line envelope is ~ 10 2 times larger and 
the emission measure is <~ 10 3 times larger than 
in the adiabatic case. 

The line cooling zone and the photoionized en- 
velope may be observed as a broad (<~ 1000 km 
s _1 ) emission line component perhaps of low in- 
tensity if compared to the narrow line caused by 
the central HII region. 

It is worth noticing that the luminosity of the 
central HII region decays rapidly after ~ 3Myr, 
when the most massive stars in a coeval cluster be- 
gin to move away from the main sequence, whereas 
the broad component conserves its luminosity be- 
ing ionized by the soft X-ray radiation and there- 
fore should be easier to detect in old (> 10 7 yr) 
objects. 

The X-ray luminosity of the star cluster wind 
is given by 

L x =4tt [ ' r 2 n 2 w A x (Z w ,T w )dr, (16) 
Jo 

where n w {r) is the atomic density distribution, 
Rx.cut is the X-ray cut-off radius where the wind 
temperature drops below 5x 10 5 K, and A X (Z, T) is 
the X-ray emissivity derived by Raymond & Smith 
(1977) in their hot-plasma code (see Strickland & 
Stevens 2000). 

5. Comparison with the observations 
5.1. The Arches cluster 

The Arches cluster is the densest and the most 
compact star cluster known in the Local Group. 
It is located within ~ 0.2pc volume at <~ 50 pc 



from the Milky Way center and contains ~ 120 
stars with masses in excess of 20 M©(Serabyn et al. 
1998). The age of the cluster is estimated within 
a range 2 - 4.5 Myr and the mass is ~ 10 4 M Q for 
an IMF with a = 1.6 and lower and upper mass 
cutoffs of lM and 100 M , respectively. 

Two sets of calculations for the Arches cluster 
wind have been presented by Raga et al. (2001) 
and Stevens & Hartwell (2003) . They differ some- 
what on the assumed input parameters. Stevens & 
Hartwell (2003) derived the total mass deposition 
rate M sc — 7.3 x 10 _4 Mq yr -1 and the average 
individual stellar wind terminal speed =2810 
km s" 1 from the Lang et al. (1999, 2001) obser- 
vations of the Arches cluster individual stars and 
adopted a Solar gas metallicity. This set of pa- 
rameters leads to a total energy deposition rate 
E sc w 1.8 x 10 39 erg s" 1 . Raga et al. (2001) as- 
sumed a lower mean individual stellar wind termi- 
nal velocity (14c = 1500 km s -1 ) and presented 
results for a star cluster with 60 identical massive 
stars, each ejecting = 10 _4 M Q yr _1 . This im- 
plies a total energy input rate E sc w 4.2 x 10 39 erg 
s _1 . Note that the adiabatic wind X-ray lumi- 
nosity between 0.3 and 8.0 keV for this set of 
parameters and Z w — 2Z Q is L x s=s 3 x 10 37 
erg s _1 , approximately two orders of magnitude 
above the value indicated by Raga et al. (2001) 
(L x r; 3 x 10 35 erg s _1 ). However we recover 
a good agreement with their results if we adopt 
= 1O~ 5 M yr -1 (instead of the cited value 
= 10~ 4 M Q y^ 1 ) for the individual stellar 
mass loss rate. Fact that indicates a missprint 
throughout their paper. 

The results of the calculations that include ra- 
diative cooling for the modified (M* = 10 M 
yr -1 ) Raga et al. (2001) star cluster model are 
presented in Figure 5a. For both sets of input pa- 
rameters (modified Raga et al., 2001 and Stevens 
& Hartwell, 2003) the Arches cluster wind evolves 
in both cases in the quasi-adiabatic regime. For 
the Raga et al. (2001) star cluster parameters the 
central temperature is approximately 3.3 x 10 7 K. 
It drops to the X-ray cut-off value (5 x 10 5 K) 
at 2.6 pc distance and to 10 4 K at 46.2pc ra- 
dius. The calculated 0.3-8.0kcV X-ray luminos- 
ity is L x = 3 x 10 35 erg s _1 , the broad emission 
line luminosities arc Ln a = 4.9 x 10 34 erg s _1 , 
L B r~( = 4.7 x 10 32 crg s _1 , respectively. For the 
Stevens & Hartwell (2003) parameters (M sc = 
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7.3 x 10" 4 M Q yr- 1 , V^a = 2810 km s^ 1 ) the cen- 
tral temperature is T c w 1.2 x 10 8 K, the X-ray 
cut-off radius is 6.9 pc and the inner boundary of 
the photoionized envelope (the 10 4 K radius) is 129 
pc. The calculated X-ray luminosity between 0.3 
and 8.0 keV and the broad emission line luminosi- 

L Hn = 6.6 x 10 34 - 



ties are L x = 10 35 erg s 1 



*erg 



s 1 and Lbt-i = 6.3 x 10 32 erg s respectively. 

5.2. The nucleus of NGC 4303 

The galaxy nucleus of the NGC 4303 belongs to 
the class of low-luminosity active galactic nuclei. 
The energy output rate from the nuclear region is 
dominated by the compact (R sc w 1.55 pc) and 
massive (M sc w 10 5 M Q ) super-star cluster. The 
Ha luminosity derived from the WHT/ISIS Ha 
flux corrected for the absorption and aperture ef- 
fects is 1.2 x 10 39 erg s" 1 (Colina et al. 2002). 
The SSC ultraviolet spectrum is best fitted by a 4 
Myr old instantaneous starburst of 10 5 M Q with 
a Salpeter IMF and 1 M Q and 100 M Q lower and 
upper stellar mass cutoffs. The thermal compo- 
nent of the unresolved-core X-ray spectrum is best 
fitted by Tw 7.5 x 10 6 K plasma with the X-ray lu- 
minosity between 0.07 keV and 2.4 keV around 

2 x 10 38 erg s" 1 (Jimcncz-Bailon ct al. 2003). 

We adopt a SSC mechanical luminosity, L sc 

3 x 10 39 erg s _1 , predicted by the Leitherer et 
al. (1999) starburst model and associate the ob- 
served hot plasma temperature with the wind cen- 
tral temperature and derive the wind terminal ve- 
locity by iterations, from the condition that the 
flow crosses the star cluster surface with the local 
sound speed. The calculated wind terminal veloc- 
ity is 715 km s _1 . Figure 5b present the temper- 
ature distribution for this strongly radiative wind 
and shows how it begins to deviate from the adia- 
batic profile (solid line) at a distance ~ 6 pc away 
from the center. It falls to the X-ray cutoff value 
at 5.9 pc and reaches 10 4 K value at 31.9 pc. The 
calculated X-ray luminosity for 0.3 and 2.0 keV 
energy range and broad emission line luminosities 
are L x = 1.3 x 10 38 erg s"\ L Ha = 1.5 x 10 36 erg 
s _1 and Lsr-y = 1-4 x 10 34 crg s _1 , respectively. 
This implies that the expected Ha broad luminos- 
ity constitutes only 0.1% of the NGC 4303 core 
Ha emission. 

Figure 5b also demonstrates how the tempera- 
ture profiles are sensitive to the adopted ejected 
gas metallicity. In the case of Z w = 3Z Q the X- 
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Fig. 5. — The stationary wind temperature dis- 
tributions, a) The Arches cluster. Solid lines are 
Stevens & Hartwell (2003) and modified Raga et 
al. (2001) solutions. The dotted lines present 
radiative solutions for the same sets of the star 
cluster parameters, b) The NGC 4303 nuclear 
super-star cluster. The adiabatic temperature dis- 
tribution is shown by the solid line. The dotted 
and dashed lines display the radiative model pre- 
dictions for solar and 3 times solar wind metal- 
licities, respectively, c) The temperature distri- 
bution for NGC 5253 supernebula critical out- 
flow. The adiabatic stationary wind solution (solid 
line) is compared with the radiative solution when 
Z w — 1.5Z Q . For larger values of Z w and the 
adopted cluster parameters the stationary solution 
does not exist. 



ray cut-off radis is 5.6 pc, the temperature drops 
to the 10 4 K value already at 14.6 pc and the calcu- 
lated wind luminosities are L x — 2.6 x 10 38 erg s -1 , 
L Ha = 1.3 x 10 36 crg s" 1 and L Bn = 1.3 x 10 34 erg 
s -1 , respectively. 

5.3. The NGC 5253 supernebula 

NGC 5253 is a nearby (3.8 Mpc) peculiar dwarf 
galaxy containing numerous massive super-star 
clusters (Meurer ct al. 1995; Gorjian 1996). An 
extraordinary compact and bright radio-infrared 
super-nebula has been discovered within the cen- 
tral star forming region of this galaxy by Turner et 
al. (2000) and Gorjian et al. (2001). The Lyman 
continuum rate required to maintain the ionization 
of the super-nebula, Nl vc = 4 x 10 52 s -1 , requires 
of a (5 — 7) x 10 5 M Q star cluster with energy depo- 
sition rate E sc « 2 x 10 40 erg s^ 1 and a Salpeter 
mass distribution having 100M Q and lM mass 
cut-off limits (see Leitherer et al. 1999). The 
mean radius of the ionizing cluster is R sc w 0.7 
pc. Using the Keck Telescope recombination line 
spectra, Turner et al. (2003) obtained a recombi- 
nation linewidth of 75 km s" 1 and concluded that 
the super-nebula gas may actually be bound by 
the gravitational pull of the super-star cluster. 

It is worth noticing that the above parameters 
imply that NGC 5253 super-star cluster is close 
to our critical energy limit (see Figure 3) if one 
assumes that the energy to mass deposition rates 
ratio is close to the {2q e /q m ) 1 / 2 ~ V^.A = 1000 
km s -1 (solid line in Figure 3). If this is the case 
a small increase in the ejected gas metallicity (to 
more than Z w = 1.5Z & ) will move star cluster 
above the critical value into the forbidden param- 
eter space. This implies that the stationary wind 
solution may not exist for this particular cluster 
and that the ejected gas is to accumulate in the 
neighborhood of the star cluster. In such a case, 
the super-nebula may consist of thermalized mat- 
ter injected by winds and supernovae, which un- 
der strong radiative cooling acquires a low sound 
speed value and is thus unable to stream at high 
velocities away from the cluster as a stationary 
wind. 

6. Conclusions 

We have developed a self-consistent stationary 
solution for spherically symmetric winds driven by 



compact star clusters taking into consideration ra- 
diative cooling. 

We have shown that stationary radiative winds 
differ strongly from their adiabatic counterparts. 
In particular we have shown that in the energy-size 
plane, there is a regime where the stationary wind 
solution is inhibited. This occurs whenever the en- 
ergy radiated away per unit volume and per unit 
time {n 2 c A(Z w ,T w )) surpasses a value of ~ 30% 
of the energy injection rate. In this catastrophic 
cooling regime, the sonic point cannot be accom- 
modated at the star cluster surface and the sta- 
tionary wind solution does not exist. Below such 
a limit the flow, despite radiative cooling, behaves 
within the star cluster volume in a quasi-adiabatic 
manner and is able to set the sonic point at the 
star cluster boundary and evolve into a stationary 
wind. 

Stationary winds driven by stellar clusters with 
an energy input rate or a size that approaches the 
critical value, establish a temperature distribution 
radically different from that predicted by the adia- 
batic solution. In these stationary wind cases the 
fast fall of temperature brings the boundaries of 
the X-ray zone, and of the line cooling zone and 
the photoionized envelope, closer to the star clus- 
ter center. This promotes the establishment of 
a compact ionized gaseous envelope which should 
be detected as a week and broad (~ 1000 km s _1 ) 
emission line component at the base of a much 
narrower line caused by the central HII region. 

Note that the threshold energy input rate ap- 
proaches an asymptotic value for large values of 
R sc (sec Figure 3). This implies that single super- 
massive star clusters are not able to generate sta- 
tionary outflows whatever their radii may be. The 
fate of the ejected material in this case remains 
unclear. A self-regulating star forming region may 
form and may keep the injected gas bound because 
of catastrophic radiative cooling or the gas may 
be blown away in a quasi-recurrent regime. The 
outflows driven by supermassive or super-compact 
star clusters should be studied with a full non- 
stationary hydrodynamic approach. 

We have speculated that the super-nebula in 
NGC 5253 seems a good example of this inhib- 
ited stationary wind regime. Radiative cooling en- 
forces a rapid drop in the sound speed value and 
the injected matter is to remain near the cluster. 
In such a case we predict that the metallicity of 
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the super-nebula is above solar, making the cluster 
lie above the threshold limit for stationary winds 
(Figures 3 and 5). 

Our calculations show that the Arches cluster 
wind seems to evolve in the quasi-adiabatic regime 
and predict the H Q and Br7 broad component lu- 
minosities around Lu a ~ 5 x 10 34 erg s" 1 and 
Lsr-y ~ 5 x 10 32 erg s _1 , respectively. 

The temperature distribution derived for the 
NGC 4303 central 1.55 pc star cluster wind is radi- 
cally different from the adiabatic temperature dis- 
tribution even for a solar wind mctallicity. The 
calculated X-ray luminosity is in reasonable agree- 
ment with the observed diffuse component lumi- 
nosity. The radiative model also predicts a com- 
pact (between 6 pc and 30 pc) broad line emission 
with Lro ~ 10 36 erg s _1 and Lsr-y ~ 10 34 erg 
s -1 , respectively. 
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thanks also to C. Munoz-Tunon for multiple sug- 
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Fe workshop. We thank prof. J. Palous for his 
comments and suggestions as a referee. We also 
appreciate the financial support given by Mexico 
(CONACYT) research grant 36132-E. 
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